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' Abstract. Recent results with the Potsdam Wolf-Rayet (PoWR) models have 

shown that Wolf-Rayet mass loss can be explained by radiative wind driving. An 
i- inspection of the galactic WR sample, however, reveals that a significant part 

D ' of the observed WR stars are in conflict with our models. This group is chiefly 

I formed by intermediate spectral subtypes. Among the population of late-type 

iy»^ . WN stars we find that the enigmatic WN 8 subtypes, which are well-known for 

' their photometric variations, are in disagreement with our models. So are there 

other driving mechanisms working in these objects? 
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\^ , 1. PoWR hydrodynamic model atmospheres 
G\ , 

I The Potsdam Wolf-Rayet (PoWR) hydrodynamic model atmospheres combine 

■ fully line-blanketed no n-LTE models with the equations of hydrodynamics (see 

^ I Grafener Sz Hamannl2 005: Ha mann Sz Grafeneil2003l : lKoesterke et all2002l : lGrafener et alJ 



In these models the wind structure {p{r), v{r), and T{r)) is computed 



in line with the full set of non-LTE populations and the radiation field in the 
O ' co-moving frame (CMF). In this numerically rather expensive approach no sim- 

^ , plifying assumptions for the computation of the radiative acceleration arad are 

^ ' made. In contrast to all previous wind models Crad is obtained by direct inte- 

gration instead of making use of the Sobolev approximation 

X ■ flrad = - / XuF^diy. (1) 

lJ ■ cJ 



In this way, complex processes like strong line overlap or the redistribution of 
radiation are automatically taken into account. The models thus describe the 
conditions in WR atmospheres in a realistic manner, in particular also at large 
optical depth. Moreover they provide synthetic spectra, i.e. they allow for a 
direct comparison with observations. 

Utilizing these m odels, we obtained the first fully self-consistent Wolf-Rayet 
wind model (Grafene r Sz Hamannll2005l l . for the case of an early-type WC star. 
Moreover we have recently exami ned the mass loss from late -type WN stars 
and its dependence on metallicity ((Grafener Sz IIamannll2006al lbll. To illustrate 
the complexity of the involved processes we show in Fig.^ the distribution of 
the radiative flux, the opacity, and their product (i.e., the radiative acceleration 
per frequency interval) within the WC model. It can be clearly seen that the 
initial FUV flux, originating from the bottom of the wind, is blocked by succes- 
sively recombining continua. This radiation is re-emitted at lower frequencies 
in the form of emission lines, forming the typical WR emission line spectrum. 
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Figure 1. Self-con s istent wind model for an early-type WC star from 
iGrafener fc Haman n ( 2005) : radiative flux (top) , mass absorption coefficient 
(middle) , and the resulting radiative acceleration per frequency interval (bot- 
tom) are plotted vs. wavelength, and wind density (as depth index). The 
wavelength ranges from X-ray to IR and the depth from the inner boundary 
at tross = 20 outward to w 1000 The opacity distribution shows succes- 
sively recombining continua due to C, O, and He in the FUV, and strongly 
overlapping line opacities in the UV. The radiative acceleration is mainly due 
to Fe-group line opacities absorbing the complex radiative flux. 
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The radiative acceleration predominantly arises from Fe-group line opacities ab- 
sorbing this radiation. These lines partly form the observable iron forest short- 
ward of 1500 A. They show a pronounced ionization structure and contribute 
by themselves to the complex 'background' flux. This coexistence of line emis- 
sion and absorption, in combination with strong overlap, is the most probable 
reason for substant ial deviations from the standard theory of line-driven winds 
( Castor et al.l IToTSl . see also Bjorkman this volume). In our WR models we 
tend to detect extr emely small or even negative force multiplier parameters (see 
iGrafener Sz Hanianni.2005i l . which can be explained by a 'statistical' line locking 
effect due to the strong overlap. 

In addition to the complex radiative processes we find a strong influence of 
wind clumping on the dynamics of O and WR star winds ( Grafener Sz HamannI 
l2003l . [2005l ). Wind clumping particularly affects the recombination rates because 
these scale with . A down-scaling of M oc 1/ y/D (where D denotes the clump- 
ing factor) thus preserves the emission measure in the atmosphere. Because WR 
atmospheres are dominated by recombination processes, the ionization structure 
as well as the emi ssion line spectrum are, app roximately, preserved under such 
a transformation ( Hamann Sz Koesterkdll998l ^. The same holds for the radia- 
tive force. The resulting radiative acceleration thus scales with ^/D under this 
transformation. 



2. Spectral analyses of galactic WR stars 

Before we discuss the results from our hydrodynamic models we give an overview 
of recent spectral analyses of galactic WR stars, based on line-blanketed PoWR 
models. In Fig. [21 we sh ow an empirica l HRD with stellar te mperatures and 
lumin osities taken from iHamann et al.l ()2006l . WN stars) and iBarniske et al.l 
( 20061 . WC stars). According to these new analyses, the galactic WR stars 
form two distinct groups in the HR diagram, which are distinguished by their 
luminosities. The first group is formed by H-rich WNL stars with luminosities 
above 10^ Lq. These stars are found to the right of the ZAMS. The second 
group are early- to intermediate subtypes with lower luminosities and hotter 
temperatures. The majority of these objects is found to be H-free. 

The observed dichotomy among the WN subtypes implies that the H-rich 
WNL stars are the descendants of very massive stars, possibly still in the pase 
of central H-burning. The earlier subtypes (including the WC stars) are less 
massive, He-burning objects. Note that distances are only known for a small part 
of the WNL sample. Some of these objects might thus have lower luminosities 
and be the direct progenitors of the earlier subtypes. 

A fundamental problem arises from the frequent detection of H-free WR 
stars of intermediate spectral subtypes, with stellar temperatures of = 50- 
100 kK. Due to their high mean molecular weight, these stars are actually ex- 
pected to lie on the He main-sequence with effective core temperatures above 
100 kK. Although the new spectral analyses with line-blanketed models tend to 
give higher values of T^, this long-standing problem still persists. Note that for 
a part of these stars the wind is so thick that the effective photosphere lies far 
above the hydrostatic radius. For these cases, T^, which is the effective temper- 
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Figure 2. Recent s pectral analyses of gal actic WR stars with l i ne-bl anketed 
models, according to iHamann et al.1 l|2006j) and lBarniske et aP l)2006i|) : sym- 
bols in light grey denote H-rich WR stars, whereas H-free objects are indicated 
in dark grey. WC stars are indicated in black. For objects with large symbols 
direct distance estimates are available ^vaji der Hucht 2001), whereas objects 
with small symbols are calibrated b y their spectral subtype. Eyolutionary 
tracks for non-rotating massive stars l|Mevnet fc Maed'e^l2003(l are shown for 
comparison. 



ature related to -R(tross = 20) only gives a lower limit for the hydrostatic core 
temperature. 



3. Hydrodynamic models for early spectral subtypes 

As mentioned earlier, our hydrodynamic model atmospheres are capable to re- 
pro duce the observed propertie s of early-type WC st ars in a self-consistent m an- 
ner (|Grafener fc Hamannll2005h . In agreement with iNugis fc Lai^I^ (|2002l ) we 
find that high critical-point temperatures (~ 200 kK) are needed to match the re- 
gion where the "hot Fe opacity peak" becomes active. According to our models, 
this also requires extremely high stellar temperatures of > 100 kK. Otherwise, 
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Figure 3. Wind acceleration within our models in units of the local grav- 
ity, vs. wind density as depth index (compare Fig.QJ. The left panel shows 
our hydrodynamic WC model with = 140 kK. Self-consistency is obtained 
throughout the complete atmosphere (flwind = vv' + g). The right panel shows 
a cooler model with = 85 kK. For this model no consistency is obtained. 
The hot Fe-peak opacities (Feix-xvii) are located so deep inside the atmo- 
sphere, that they are missing in the wind acceleration region. 



as illustrated in Fig.Ol the hot Fe-peak moves so deep inside the atmosphere that 
no line opacities are available for the acceleration of the wind base. For inter- 
mediate spectral subtypes with = 70 -100 kK our models thus predict no 
WR-type mass loss. The existence of H-free WR stars in this temperature range 
is thus ruled out by the theory of stellar structure and our wind models. 

Nevertheless, we have seen that such objects are commonly observed (see 
Fig. 12). The additional extension of their envelopes, as well as their strong mass 
loss, can not be explained by standard assumptions. We thus conclude that 
other processes might be important for these objects. In particular, this means 
that their stellar winds might not be purely radiatively driven. According to our 
models the lack of wind acceleration occurs at the bottom of the stellar wind. 
Fast rotation seems to be a good candidate for an additional support of the wind 
in this region. 



4. Hydrodynamic models for late spectral subtypes 

In a recent work we have investigated the prop erties of luminous, H-rich WNL 
stars with hydrodynamic PoWR models (Grafener & Hamann 200(^. The 
most important conclusion from that work is that WR-type mass loss is pri- 
marily triggered by high L/M ratios or, equivalently, Eddington factors Fe = 
XeLi^/ATTcGM-i, approaching unity. Due to the increase of the density scale height 
at low effective gravities, the formation of optically thick winds seems to be 
strongly favored for such objects. Note that high L/M ratios are expected for 
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Figure 4. Wind models for galactic WNL stars: w ind densities compared to 
observed values according to lHamann et all i)200(t() . In the i^t-T* plane, our 
models (solid curves) reproduce the upper part of the observed WNL sample 
(crosses), corresponding to WNL stars with low wind densities (see text for 
the definition of Rt). The objects belonging to the enigmatic class of WN8 
spectral subtypes are indicated in grey. 



very massive stars and for He-burning objects, giving a natural explanation for 
the occurrence of the WR-phenomenon. 

In Fig. 1^ we show a quantitative comparison of the results from our hydro- 
dynamic models wit h observations of the galactic WNL sample (according to 
iHamann et ani2nn^ l. For this comparison we employ the transformed radius 



Rt — Ri, 



DM 



2500kms-i / lO-^Mgyr- 



2/3 



(2) 



which is a luminosity-independent measure of the wind density. Models with 
the same -Rt are scale invariant, and show a very similar spectral appearance for 
given and given surface chemistry (see Hamann & : Koesterke 1998). In the 
Rt-T-k plane, our models with Eddington factors Fe = 0.55 and 0.67 reproduce 
the observed WNL stars with large i?f These are the WNL stars with relatively 
low wind densities. The resultant mass loss rates depend strongly on the stel- 
lar temperature (with M oc T^'^'^), nicely reflecting the observed sequence of 
spectral subtypes from WN 6 to WN 9. 

The high wind density objects, however, are not reproduced by our models. 
Notably, these objects are dominated by the enigmatic WN8 subtypes, which 
tend to show strong e mission lines in com bination with photometric variations 
(e.g. iMarchenko et al.il998. : ■Lefevre et al.ir2005 ) . Because of the apparent depen- 
dence of M on Fe in Fig.|^ one might think that higher mass loss rates can easily 
be obtained by a further increase of Fg. In our models, however, this leads to a 
point where the radiative acceleration exceeds gravity in the hydrostatic layers. 
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For the relatively cool WNL stars this situation occurs due to the hot Fe- 
peak which causes an inward increase of the Rosseland mean opacity from a 
certain point on. If gravity is exceeded by the resulting radiative force, a sta- 
tionary solution can only be obtained if the excess acceleration is compensated 
by an outward increase of the gas pressure, i.e. a density inve rsion. I i i real - 
ity this will presumably lead to a non-stationary situation. .Dorfi et al.l ( 20061 ) 
have recently shown that the observed variability in WN 8 subtypes can indeed 
be explained by instabilities which are caused by the hot Fe-peak. The ob- 
tained timescales aro und 0.5 d are in rou gh agreement with the observed period 
of 9.8 h for WR123 ((Lefevre et al.ll2'005^ . and with the dynamical timescale of 
the wind (-R^/voo ~ 4h). We thus conclude that the observed pulsations will 
most probably influence the mass loss of these objects, i.e., their winds are not 
purely radiatio n driven. Note that the non-detection of X-ray emission from 
WN 8 subtypes (jOskinovall2005h might also be a hint on a different wind driving 
mechanism for these objects. 

Another reason for our difficulties with the WN 8 subtypes might be related 
to the fact that the distances of these objects are only poorly known. They might 
thus fall into the category of core He-burning post-RSG stars with much lower 
luminosities. A detailed investigation of the influence of on our wind models, 
however, remains a topic of future investigations. 

Finally, due to their proximity to the Eddington limit, our models allow to 
estimate the masses of WNL stars. Because of the strong dependence of M on the 
L/M ratio, we can give reliable mass estimates, that only depend on the adopted 
distance to the object. For the H-rich, weak-lined WNL stars our models tend 
to give very high masses, in agreement with very massive stars in a late stage of 
the H-burning phase. For the case of WR22, a weak-lined WN 7h subtype with 
= 10^ '^ -Lq, we derive a s t ellar mass of 78 Mq. This value is in agreement 
with the work o f lRauw et aP (|l99fil 'l. who derived 

sin i = (72 it 3) Mq from 

orbital elements. 



5. Conclusions 

Our hydrodynamic atmosphere calculations for WR stars have shown that WR- 
type stellar winds can be driven alone by radiative acceleration. In particular 
we found that the proximity to the Eddington limit is the primary trigger of the 
enhanced mass loss of WR stars. Nevertheless, difficulties remain with H-free 
objects of intermediate spectral subtype, and late-type WN stars with strong 
stellar winds (mostly of spectral type WN 8). In both cases, problems are caused 
by the gap between the 'hot' and the 'cool' Fe opacity-peaks. For the intermedi- 
ate subtypes this gap causes a deficiency of line opacities at the wind base. This 
problem could be overcome by an additional support due to fast stellar rotation. 
For the WN 8 stars, the inward increase of the radiative force due to the hot Fe 
opacity-peak causes problems in the hydrostatic layers, presumably leading to 
the observed variability. These oscillations most probably affect the mass loss 
of these stars. 
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